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ABSTRACT 

We present ISOPHOT observations between 60 and 200 /ini and a near-infrared ex- 
tinction map of the small intermediate-density cloud LDN 1780 (Galactic coordinates 
1=359° and b = 36.8°). For an angular resolution of 4', the visual extinction maximum 
is Av = 4.4 mag. We have used the ISOPHOT data together with the 25, 60 and 100 
/im IRIS maps (Miville-Deschenes & Lagache 2005) to disentangle the warm and cold 
components of large dust grains that are observed in translucent clouds (Cambresy 
et al. 2001, del Burgo et al. 2003) and dense clouds (del Burgo & Laureijs 2005). 
The warm and cold components in LDN 1780 have different properties (temperature, 
emissivity) and spatial distributions, with the warm component surrounding the cold 
component. The warm component is mainly in the illuminated side of the cloud facing 
the Galactic plane and the Scorpius-Centaurus OB association, as in the case of the HI 
excess emission (Mattila & Sandell 1979). The cold component is associated with the 
^•^CO (J=l-0) Hne integrated {W13), which trace molecular gas at densities of ~10^ 
cm^^. The warm component has a uniform colour temperature of 25±1K (assuming 
P = 2), and the colour temperature of the cold component slightly varies between 
15.8 and 17.3 K {(3 = 2, AT=0.5K). The ratio between the emission at 200 fim of 
the cold component (/^(200)) and Ay is /^(200)/Av/=12.1±0.7 MJy sr'^ mag-^ and 
the average ratio T2oo/Av = (2.0 ± 0.2) x 10~^ mag^^. The far infrared emissivity of 
the warm component is significantly lower than that of the cold component. The Ha 
emission (/^(Hq!)) and Ay correlate very well; a ratio I^{lla)/Av = 2.2±0.1 Rayleigh 
mag~^ is observed. This correlation is observed for a relatively large range of column 
densities and indicates the presence of a source of ionisation that can penetrate deep 
into the cloud (reaching zones with optical extinctions Ay of 2 mag). Based on mod- 
elling predictions we reject out a shock front as precursor of the observed Ha surface 
brightness although that process could be responsible of the formation of LDN 1780. 
Using the ratio Iij{Ila)/Av we have estimated a ionisation rate for LDN 1780 that 
results to be ~10^^® 7 s^^. We interpret this relatively high value as due to an en- 
hanced cosmic ray flux of ~10 times the standard value. This is the first time such an 
enhancement is observed in a moderately dense molecular cloud. The enhancement in 
the ionisation rate could be explained as result of a confinement of low energy ('--^lOO 
MeV) cosmic rays by self generated MHD waves (Padoan & Scalo 2005). The origin 
of the cosmic rays could be from supernovae in the Scorpio-Centaurus OB association 
and/or the runaway C Ophiuchus. The observed low ^^CO abundance and relatively 
high temperatures of the dust in LDN 1780 support the existence of a heating source 
that can come in through the denser regions of the cloud. 

Key words: ISM: clouds, dust, extinction - infrared: ISM. 



1 INTRODUCTION 
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The Lynds Dark Nebula LDN 1778/1780 (Lynds 1962) is a 
moderately dense (translucent) region with Galactic coor- 
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dinates 1 = 359° and b=36°.8. This compact nebula, here- 
inafter only referred as LDN 1780, is located at h « h© (=15 
pc) + 66 pc (assuming a distance of llOilO pc, Franco 1989) 
from the Galactic midplane. It is 1.28 pc x 0.96 pc in size, 
slightly elongated in the E-W direction on the sky. It shows 
a well defined boundary towards NW and is more diffuse on 
the SE side, which faces the Galactic Plane. 

LDN 1780 was studied together with other nearby 
clouds (LDN 134, LDN 183, LDN 169, MBM 38 and MBM 
39) by Laurcijs ct al. (1995, hereinafter LFHMIC95) from 
carbon monoxide observations at 2.7 mm of CO(J=1-0), 
^^CO(J=1-0) and ^^CO(J=1-0)S 12, 25, 60 and 100 /xm 
/i?^S'data, and blue extinction {Ab) from star counts. They 
found that the ^^CO(J=1-0) abundance in LDN 1780 is fa4 
times lower than in the opaque clouds LDN 134, LDN 183 
and LDN 169. The 100 jim. emission excess A7i.(100) = 
7^(100) - /^(60)/e (where 7^(60) and 7^(100) are respec- 
tively the emissions at 60 and 100 jira and O is the average 
ratio 7i,(100)/7^(60), which is 0.27 for LDN 1780 and 0.21 
for the rest of the clouds) corresponds to the emission from 
the "classical" big grains (without the contribution from the 
very small grains). LFHMIC95 found that there is a good 
correlation of A7,,(100) and the ^^CO(J=1-0) line integrated 
(Wis). Although the ratio Ah{lQQ)/Wv.i in LDN 1780 is 
significantly higher than in the other clouds, it presents a 
similar ratio A7„(100)/Ab. LFHMIC95 concluded that the 
anisotropic UV radiation field in the complex affects the 
density distribution, especially on the illuminated side of 
the clouds with a bigger halo, but also in the densest, well- 
shielded regions, shifted with respect to the cloud centres to 
the opposite side of the illuminated halos. Only LDN 1780 
lacks a diffuse halo and the radiation field is still affecting 
the molecular content given the low opacity of the cloud 
(LFHMIC95). 

LDN 1780 is near Radio Loop I (Berhuijsen 1971), likely 

just at the boundary of the Local Bubble (a Galactic region 
of radii 65-150 pc around the Sun of low density and hot gas 
with a temperature of 10® K; see Lallemcnt ct al. 2003 and 
references therein) . Radio Loop I is thought to be the limb of 
the neighbouring hot bubble Loop I, a superbubble centred 
on the Sco-Ccn OB association (distance of 120-160 pc, de 
Geus, de Zeeuw & Lub 1989) that is expanding at 10-20 
km s~^ due to stellar winds and supernovae (de Geus 1992). 
X-ray observations at i, | and 1.5 KeV support that LDN 
1780 is located at the boundary of the Local Bubble and 
show that the cloud is not affected by any possible hot gas 
flowing from the hotter bubble Loop I to the Local Bubble 
(Kuntz, Snowden & Verier, 1997). 

LDN 1780 is the only translucent cloud in which has 
been observed Extended Red Emission (ERE) (Chlewicki & 
Laurcijs 1987) so far. The ERE in this cloud has a broad, 
unstructured band, with a peak near 7000 A. This peak's 
wavelength is longer than expected, and it could be due to 
silicon nanoparticles in a environment with relatively high 
electron density and lower UV photon density in comparison 
with the diffuse interstellar medium (Smith & Witt 2002). 

Observations in the 21 cm HI line and the 18 cm OH 
lines of LDN 1780 were carried out using the Effelsberg 100 
m telescope by Mattila & Sandell (1978). They concluded 



1 No detection of 18CO(J=1-0) in LDN 1780 was achieved. 



that the atomic and molecular gas arc spatially connected 
with cacii other and the dust cloud. Observations of LDN 
1780 at 9 cm CH (Mattila 1986), 6 cm H2CO (Clark & John- 
son 1981) have been performed too. Molecular and HI excess 
emission spectra show a peak intensity at the local-standard- 
of-rest velocity vlsr = 3.5 km s~^ and broad blue shifted 
profiles. Based on CO(J=1-0), "CO(J=1-0) and C^®0(J=1- 
0) observations Toth et al. (1995)^ (hereinafter THLM95) 
claimed the observed cometary like structure and morphol- 
ogy of LDN 1780 is result of shock fronts from the Sco-Cen 
OB association, as well as an asymmetric UV radiation from 
the Galactic disk and the Upper Scorpius group. 

LDN 1780 presents^ a unique "CO(.1=1-0) core with a 
regular morphology (LFHMIC95, THLM95), a narrow ve- 
locity component of wO.6 km s"^ (LFHMIC95; 0.52 km 
s"^ from THLM95), an average density of 10^ cm"^ (LFH- 
MIC95; a lower value of 0.6 10^ cm~^ is obtained from its 
column density and size determined by THLM95), and a to- 
tal mass of almost 10 Mq (LFHMIC95; 8.3 Mq according 
to THLM95). It displays a r"^ density distribution (LFH- 
MIC95) and it is in virial equiUbrium (THLM95). The total 
meiss in the whole ^^CO(J=1-0) emitting region is 18 M© 
(THLM95), and the estimated mass of the cloud envelope is 
1.8-3.6 Mq (THLM95). The total mass of LDN 1780 is then 
«20 Mq. 

In this paper we present observations of LDN 1780 ob- 
tained using ISOPHOT in the spectral range 60-200 yitm 
that trace the emission from very small grains (VSGs) and 
big grains (BGs). These far- infrared data have been com- 
plemented with Two Micron All Sky Survey (2MASS) near- 
infrared photometry (Cutri et al. 2003), the 25, 60 and 
100 /um IRIS maps (Miville-Deschenes & Lagache 2005), 
Hq maps (Finkbeiner 2003), and CO molecular line mea- 
surements (LFHMIC95). We have studied the distribution 
and the properties (temperature and emissivity) of the dust 
grains in the cloud, and compared some of these properties 
with the molecular gas emission, the Ha emission and their 
distributions. An extinction map is derived from 2MASS 
photometry, and converted to optical extinction Ay assum- 
ing a standard extinction law. This map is used as a tracer of 
the column density of dust at short wavelengths. The warm 
and cold components of large grains (see del Burgo et al. 
2003, hereinafter Paper I) are separated following a similar 
method to that used by del Burgo & Laureijs (2005) (here- 
inafter Paper II). The extinction map is compared with the 
emission and the optical depth at 200 /im of the cold com- 
ponent, to investigate emissivity variations of the big dust 
grains across the cloud, as carried out in previous studies 
(Paperl; Paper II). The Ha surface brightness in the cloud 
is corrected from extinction, and a cosmic ionisation rate is 
derived from the ratio 7,y(Ha)/^v. We discuss the nature 
of the ionisation source to explain the presence of Ha emis- 
sion in moderately extinguished regions {Ay ~2 mag) in the 
cloud. We have also studied the correlations of the far in- 
frared dust tracers with molecular gas measurements of CO 
at 2.7 mm. The results for this intermediate density cloud 



^ A faint C^®O(J=l-0) emission spectra were obtained in few 
positions only. 

^ A distance of 110 pc for LDN 1780 was assumed to determine 
the following parameters. 
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are compared with other studies of translucent and dense 
regions. Section 2 describes the data processing. Sections 3 
and 4 are respectively devoted to the properties of the two 
big dust grain components and the ionisation process in the 
cloud. Finally, Section 5 presents the conclusions. 

2 ACQUISITION, PROCESSING AND 
PREPARATION OF THE DATA 

2.1 ISOPHOT observations 

A nearly squared 38' x 38' (equivalent to 1.2 pc xl.2 pc 
for a distance of 110 pc) region centred at the equatorial 
coordinates 02000=15'' 40'"09.5" and <52ooo=-7°12'28.7" and 
a stripe-shaped region centred at 02000 «15''39'"50" and 
feooo «-7°ll', were mapped with ISOPHOT (Lemke et al. 

1996) , an instrument on board of ESA's Infrared Space Ob- 
servatory {ISO, Kessler et al. 1996). The stripe- like region 
(~38 arcmin along E-W; ~3 arcmin along N-S) was observed 
using the Astronomical Observation Template PHT22 in 
raster mapping mode with the array detectors ClOO (3x3 
pixels; 46" pixeP^) and the 60, 70 and 100 ^m filterbands, 
and also with C200 (2x2 pixels; 92" pixeP^) and the 120, 
150 and 200 fj,m filterbands. The squared region was ob- 
served in the same observing mode with ClOO at 100 
and with C200 at 200 ^m, and it nearly covers LDN 1780. 
For an overview of the different observing modes see the 
ISOPHOT Handbook (Laureijs et al. 2003). 

Table lists the log of the observations, with the ISO 
identification (ISOid), the equatorial coordinates (right as- 
cension and declination J2000.0) of the raster centre, the 
reference wavelength of the ISOPHOT filter band (Are/), 
the raster dimensions (NxM), the distances between adja- 
cent points and lines (dN and dM, respectively), and the 
exposure time (texp) per sky position. The redundancy is 
generally 2 or 1. There are also some small gaps that are 
regularly interleaved along the stripe-like region when sur- 
veyed with the ClOO array detector. The total exposure time 
per pixel is obtained multiplying the redundancy and texp- 
In addition to the observations of the object, two measure- 
ments of the internal Fine Calibration Source (FCS) were 
obtained with the same filter, detector and aperture just 
before and after each map. 

The data reduction was performed with the ISOPHOT 
Interactive Analysis software PIA VIO.O (Gabriel et al. 

1997) . Signals corresponding to the same sky position were 
averaged after correcting for glitches. The slow response or 
transient behaviour of detector C200 for fiux changes is gen- 
erally visible only on the first few raster steps of the maps, 
but it did not affect the main features of these maps. Mem- 
ory effects due to crossing bright sources are negligible. All 
standard signal correction steps (reset interval correction, 
dark signal subtraction, signal linearisation) were applied. 
For flux calibration the detector's actual response was deter- 
mined from the second FCS measurement, which is observed 
just after the object. Unless explicitly stated no colour cor- 
rection was applied on the data. 

We do not correct from the by-passing sky light (del 
Burgo, Heraudeau, Abraham 2003). Note, however, our 
analysis minimises the impact of any undesirable systematic 
offset, for instance, using pixel-pixel correlations to derive 
average emission ratios (see Paper I). 



The first quartile normalisation flat-fielding method of 
PIA was used to correct for the remaining responsivity dif- 
ferences of the individual detector pixels. 

Fig.[Tland|lshow the ISOPHOT emission maps of LDN 
1780 for the stripe-like and the square regions, respectively. 
Maps in Fig. |21 were smoothed to the angular resolution of 
200 ISOPHOT emission map (see Paper I). These maps 
have been zero-level calibrated following the procedure de- 
scribed in Sect. 2.4. 



2.2 Near-infrared extinction 

We built the extinction map of LDN 1780 following the 
method described in Cambresy et al. (2002) using the 
2MASS H - Ks colour excess of 3400 stars in a 80' x 80' 
field. The extinction is obtained following the expression: 

A.= [^-^Y\iH-K.)+Z... (1) 

where {Ah/ Ay - AkJAv)~^ = 15.87 (Rieke & Lebof- 
sky 1985) and Zco\ ~ —1.9 mag. The zero point for the 
extinction calibration, Zcoi, corresponds to an intrinsic star 
colour of H — Ks = 0.12 mag, measured outside the cloud. 
Any diffuse extinction in the field is therefore removed by 
construction. The extinction is obtained by taking the me- 
dian value within adaptive cells that contain a constant 
number of 3 stars. This produces a map with a varying reso- 
lution, increasing with the star density. The resulting map is 
convolved with an adaptive kernel to obtain the final maps 
at the uniform resolutions of 4' and 6' (see Fig.EJ. Although 
one can naively think this approach is equivalent to directly 
use an uniform grid of 4' or 6', it is fundamentally different. 
Because the ISM has small-scale structures, variations of the 
extinction within a cell are expected. The direct consequence 
of the sub-cell structures is that stars are not uniformly dis- 
tributed within a cell, but are preferentially detected in the 
lower extinction region of a cell. In the classical method with 
large cells, the few stars detected at higher extinction do not 
contribute significantly to the colour median of a cell. The 
weight of the low-extinction part of a cell is increased and the 
extinction measurements are biased towards small values. In 
the adaptive method with only 3 stars, cells have the mini- 
mum size and the final 4 or 6' resolution is obtained through 
a convolution. The weight of the low extinction region does 
not dominate anymore since these cells are smaller than at 
large extinction. Using adaptive cells and a convolution give 
a weight to the individual reddening measurements that is 
proportional to the surface of the voronoi cell defined by 3 
stars (the voronoi cell for a single star represents all points 
equidistant of the central star to the surrounding objects). 
This technique integrates both the star colour and density, 
which optimises the extinction estimate by reducing the bias 
due to the inhomogeneous star distribution within a cell (see 
Cambresy et al. 2005 for more details). The statistical un- 
certainties in the extinction maps are 0.4 mag and 0.5 mag 
for the 6' and 4' resolution, respectively. At these angular 
resolutions, the 6' and 4' extinction maps respectively reach 
maxima of 3.9 mag and 4.4 mag in the central region of LDN 
1780. 
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Figure 1. From top to bottom, from left to right: Emission at 60, 70, 100, 120, 150 and 200 ^im smoothed to the angular resolution 
the 200 /im emission corresponding to the stripe region, which is overplotted on the 200 /um emission map corresponding to the square 
region. / and S labels respectively stands for the illuminated and shadowed regions separated by dashed lines from the rest of the stripe 
(see Sect. 2.4.1). 
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Table 1. ISOPHOT observations of LDN 1780. 



aj2ooo <5 ,72000 ISOi^ Arcf NxM dN, dM Area texp Remark 

[hh mm ss] [° ' "] [^tm] [arcsec, arcsec] [arcmin-^] [s] 



15 39 40.93 


-7 10 35.04 


43100206 


60 


16x2 


180, 90 


206 


34 


stripe 


15 39 40.93 


-7 10 35.02 


43100207 


70 


16x2 


180, 90 


206 


34 


stripe 


15 39 40.93 


-7 10 35.01 


43100208 


100 


16x2 


180, 90 


206 


34 


stripe 


15 39 58.55 


-7 11 39.14 


43100209 


120 


13x2 


180, 90 


249 


34 


stripe 


15 39 58.55 


-7 11 39.12 


43100212 


150 


13x2 


180, 90 


249 


34 


stripe 


15 39 58.55 


-7 11 39.11 


43100210 


200 


13x2 


180, 90 


249 


34 


stripe 


15 40 09.50 


-7 12 28.74 


43100630 


100 


18x24 


135, 90 


1495 


16 


square 


15 40 09.51 


-7 12 28.69 


43100629 


200 


13x13 


180, 180 


1469 


17 


square 




-10 10 -10 10 

Ax [orcmin] Ax [orcmln] 

Figure 2. Left: ISOPHOT 100 fim emission smootlicd to tlic angular resolution of the ISOPHOT 200 fim emission; Right: ISOPHOT 
200 /im emission. Note the thin top-left arrow points to the Centaurus-Scorpius OB association. The NE axes are displayed on the 
bottom-right corner (arrow points to North). 



2.3 Additional data 

In order to complement the ISOPHOT dataset of LDN 1780, 
the following additional data were included: 



(i) IRIS maps (Miville-Deschencs & Lagache 2005) at 25, 
60 and 100 /im, with angular resolutions (noise levels in MJy 
sr"^) of 3.8' (0.05), 4.0' (0.03) and 4.3' (0.06), respectively. 

(ii) Carbon monoxide observations at 2.7 mm, (CO and 
"CO) from LFHMIC95, with an angular resolution of 2.7'. 
In particular we have used the velocity-integrated (line in- 
tegrated) emission maps of CO (W12) and ^^CO (W13). 

(iii) Ha surface brightness map, with an angular res- 
olution of 6', belonging to the H-alpha Full Sky Map 
(Finkbeiner 2003). The absolute calibration of the compos- 
ite Ha emission map is based on the undersampled WHAM 
survey (1 deg spacing; 1 deg diameter top-hat beam). The 
error in a 40'x40' Ha emission map of LDN 1780 is ~0.4 
Rayleigh (1 Rayleigh=^ 7 cm~^ s~^ sr~^). 



2.4 Data preparation 

We have applied different analysis approaches for the data 
analysis of the stripe-like region and squared field data. In 
the remainder of this paper, wo shall use the symbol Iv{\) 
(in MJy sr~^) to refer to the surface brightness, where A=25, 
60, 70 and 100, 120 and 200 fim corresponds to the refer- 
ence wavelengths of the IRIS and ISOPHOT filterbands. 
A distinction between IRIS and ISOPHOT data will be no- 
ticed if necessary. /^(Ho) refers to the Ha surface brightness. 
A/i,(A) refers to the error at the reference wavelength A. 

2.4-1 Far-infrared emission correlations: stripe-like region 

Following the analysis of Paper I, we have derived colour 
variations from the analysis of the correlations of Iv{X) 
(A <200 /im) with /,y(200) corresponding to the stripe-like 
region. As described in Paper I, any colour variation will 
result in a deviation from linearity in the correlations. The 
method offers some advantages (see Paper I for a detailed 
discussion). Slopes of correlations are direct measurements 
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Figure 3. Extinction map and contours in Ay units. Pixel size 
of 4' X 4'. Coordinates [0,0] correspond to a2000=15''40™11.6'' 
and ,52000=- 7°12'13.5". 

of the colour ratios and are independent of any emission 
offset; we therefore did not apply any zero- level correction 
to the stripe-like region. In addition, slopes are much more 
accurate than the absolute surface brightness, removing a 
significant amount of the random noise of individual pix- 
els. In particular, when temperature variations are within 
the uncertainties for a certain region, the use all indepen- 
dent points improves the accuracy and the method applied 
is optimal. 

All the stripe-like region data were obtained using 
ISOPHOT. Data corresponding to the ClOO camera were 
smoothed to the angular resolution of the 2x2 pixel C200 
camera. Fig. |1| shows the correlations for the emissions at 
60, 70 and 100, 120 versus emission at 200 ^m. Correlations 
for two separated regions of LDN 1780 have been consid- 
ered (see Fig.0 bottom): a region in the shadow part of the 
cloud and a similar-in-size region in the illuminated side of 
the cloud, which faces the Galactic plane and the OB as- 
sociation. The instrumental noises of these emissions have 
been determined from the pixel to pixel response variations 
derived from the PCS measurements (see Table |5J. 

2.4- S Map homogenisation: squared field 

The maps of IRIS, ISOPHO T, the carbon monoxide lines as 
well as the surface brightness map for the squared field 
have different angular resolutions. For a proper comparison 
the maps were smoothed to the worst angular resolution''. 

* We excluded the HI excess emission map for this comparison. 



For the analysis of far-infrared, extinction and molecular 
data, maps were smoothed to the angular resolution of the 
100 fim IRIS map of 4.3' (corresponding to 0.14 pc for a 
distance of 110 pc). For a comparison with the Ha sur- 
face brightness with FWHM«6', the maps were smoothed 
to mimic that angular resolution. 

In Fig. the 25, 60 and 100 ^lm IRIS maps (see Fig. 
top, middle-left) and 200 ^m ISOPHOT m&p (see Fig.^l 
middle-right) at the angular resolution of 4.3' are shown. 

The median instrumental noise in 7,^(25), Ii,{60), 
1^(100) reported by Miville-Deschenes & Lagache (2005) is 
0.05, 0.03 and 0.06 MJy sr~\ respectively. For ISOPHOT, 
the pixel to pixel response variations derived from the FCS 
measurements (see Sect. 2.1) have been used to determine 
the instrumental noise. Table |5| shows the values of the in- 
strumental noise at the different wavelengths for the obser- 
vations. Data were smoothed to the 4.3' angular resolution. 
When data from different instruments are combined we have 
also considered the relative errors between different filter- 
bands. In Paper II we used a 5% between ISOPHOT filter- 
bands and a pessimistic 30% for the relative error between 
IRIS 100 fim and the ISOPHOT 200 fim. Here, we have cor- 
related the IRIS and smoothed 4.3' ISOPHOT 100 fim maps 
(see Fig.|^ bottom-left) to quantify this relative error in LDN 
1780. We found an excellent correlation (Pearson correlation 
coefficient of PCC=0.97) between these maps and a relative 
error of only 1.6%. This value has been used as reference 
to derive the relative errors between IRIS and ISOPHOT 
filterbands assuming a 5% between different filterbands of 
the same instrument. 



2.4-3 Separation of the warm and cold components 

The emission maps of the squared field in LDN 1780 at 25, 
60, 100 and 200 /im show distinct morphologies that likely 
indicate variations in the abundance distribution across the 
cloud. Note, for instance, the different positions of the emis- 
sions peaks at different wavelengths in Fig. |3 See also Fig. 
0for intermediate wavelengths. 

In order to distinguish between the different large dust 
grain components in LDN 1780 we have performed a first 
attempt to obtain the emission of the cold component at 
100 j-im (see Laureijs, Clark & Prusti 1991; Abergel et al. 
1994), and at 200 ^m (see Paper II) using the equation: 

AL{X){x,y) = L{\)(x,y) - Lm{x,y)/Q{60,X), (2) 

where 0(60, A) is the average colour 1^(60)//,/ (A) of 
LDN 1780. We have determined 6(60, 100)=0.27 from the 
IRIS maps, which is the same value found by LFHMIC95 
from previously realised IRAS data. Our value of 0(60, A) 
was determined using values of 1^(100) below 90% of the 
peak of 1^(100). We found ©(60, 200)=0.13 using values of 
Ii,{200) below 90% of the maximum. The correlation be- 
tween A/^(100) and 1^(200) has a PCC of 0.7. 

The previous approach has been traditionally used to 
separate the BGs that are associated to the molecular re- 
gions and are free from VSG emission from the warm diffuse 
ISM. Although they were originally called warm and cold, 
these components separate two grain populations rather 
than two temperatures. In addition recent results obtained 
at longer wavelengths using far-infrared and submillimeter 
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Figure 4. Surface brightness /^(A) versus /,y(200) for wavelengths of 60, 70, 100, 120 and 150 fim. White and blaek circles respectively 
correspond to the illuminated region and shadowed region of LDN 1780. Straight lines that fit the data points are overplotted (see Sects. 
2.4.1 and 3.1.1). 
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Figure 5. Top-left: IRIS 25 /um emission map; Top-right: IRIS 60 /im emission map; Middle-left: IRIS 100 ;tiri ciuissioii map: Middle- 
right: ISOPHOT 200 jxm emission map. All maps with angular resolution of 4.3'. The NE axes (arrow points to North) are displayed 
on the bottom-right corner of each map. The thin and thick lines respectively correspond to the 1:1 and observed relations; Bottom-left: 
ISOPHOT 100 iMa emission versus IRIS 100 /im emission. Note emission are not colour-corrected. The thin and thick lines respectively 
correspond to the 1:1 and observed relations; Bottom-right: '^^CO(J=1-0) line integrated. Units of K Km s~^. 
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Table 2. ISOPHOT instrumental noise. 



A/^(60) 
[MJy sr-l] 


A/„(70) 

[MJy sr-l] 


A/„(100) 
[MJy sr-i] 


A/^(150) 
[MJy sr~i] 


A/„(200) 
[MJy sr~i] 


Remark 


0.4 


0.35 


0.3 


0.55 


0.5 


stripe at 1.5' 






0.45 




0.76 


square at 1.5' 






0.16 




0.27 


square at 4.3' 



instruments indicate the presence of two temperature com- 
ponents of BGs, with temperatures of «16-22 K and ~12 K, 
respectively, and different emissivities (Cambresy et al. 2001, 
Papers I and II). The designation of warm and cold compo- 
nents is still used by the community to refer to this separa- 
tion. 

A second approach has been used to carry out the dust 
component separation. Before introducing this approach, we 
give a brief scientific justification for its use, which will be 
discussed in more detail in Sect. 3.6. Assuming a single EG 
component, the observed mean ratio /^( 100)//^ (200) would 
indicate a temperature of ~17K in the cloud, which is rela- 
tively high compared with other translucent clouds (see Pa- 
per I). Also, this temperature is significantly lower than the 
temperature of ~13K derived from the A7i,(100)//„(200) 
ratio using a modified black body with a power-law emissiv- 
ity p—2. This discrepancy in temperature can be interpreted 
as due to the presence of two EG component instead of a 
single one. This could explain other findings in the cloud. It 
is observed that the 25 /xm emission resembles properties of 
both the 12 and 60 fim emission maps. In addition, the av- 
erage colour 0(60, 100)=0.27 is relatively high with respect 
to the standard value in the ISM of 6(60, 100)=0.21. LFH- 
MIC95 suggested that the high ratio 9(60, 100) in LDN 1780 
might be due to a high contribution of VSGs, but could be 
alternatively interpreted as due to the presence of a warm 
component of EGs. The enhanced stellar radiation field in 
LDN 1780 (LFHMIC95, THML95) due to the Galactic plane 
and a close OB association could produce that the warm and 
cold components present relatively high temperatures com- 
pared with those found in Cambresy et al. (2001), Papers 
I and II. A significant contribution to the emission at 60 
Hm from warm BGs is therefore expected. The 25 fim emis- 
sion, which presents a high S/N ratio, becomes the suitable 
template of the VSG emission. 

Hereinafter, when referring to the emission of the warm 
and cold components we respectively use the superscripts w 
and c. We determined the emission of the cold component 
at an intermediate wavelength (60, 100 and 200 /im) from 
the following equation: 

i:{X)ix,y) = I4X)(x,y) - c(\,25) X L(25){x,y), (3) 

where c(A, 25) is a constant value derived from the cor- 
relation J^(A)-/„(25). It corresponds to the slope of this cor- 
relation, i.e. /i,(A)//^(25), for values of /^(A) below the 90% 
of its maximum. This is performed to reject out possible 
outliers in the emission maxima surroundings. 

In order to determine the emission of the warm compo- 
nent at wavelengths A =25, 60 and 100 nm we applied the 
equation: 

CiX)ix,y) = L(\){x,y) ~ w{X,200)xi:i200){x,y), (4) 



where «;(A,200) is the ratio /„(A)//S(200). We have 
simply used /i/(200) instead of 1^(200) when determining 
/™(25). 

Note this second approach provides emission maps at 
wavelengths 60, 100 and 200 fim for the cold component. 
The 100 and 200 fj,Tn emission maps will be used to derive 
colour temperature and optical depth maps in the Sections 
3.2.2 and 3.3, respectively. For the warm component emis- 
sion maps at 25, 60 and 100 fim are obtained. The latter two 
maps have been used to determine the colour temperature 
and optical depth of the warm component (see sections 3.2.3 
and 3.3, respectively). 

2.4-4 Zero level calibration 

We have fixed to zero the minimum value of /^(lOO) in the 
field of view of ISOPHOT. The region is representative of 
the local background for LDN 1780. Then, the pixel-pixel 
correlations of 1^(60) with /^(lOO) and /^(200) with /^(lOO) 
where used to derive the offsets of the emissions at 60 and 
100 fim with respect the emission at 200 ^m. We have es- 
timated an error of A/^(200)=l MJy sr"^ due to the zero 
level uncertainty. 

For the warm component we have fixed the minimum of 
7™ (100) to zero. Then, the zero levels of the emissions at 25, 
and 60 fim and their corresponding errors were determined 
from the pixel-pixel correlations with the corrected 1^(100). 
We have computed an error at the reference wavelength of 
100 ^lm of AC(100)=0.15 MJy sr"^ due to the zero level 
uncertainty. 

We were consistent with these definitions for the zero 
level calibration when applying the method described in the 
previous section, i.e, when removing the emissions at 25 and 
200 fj,m emissions. We used the zero-level errors together 
with the relative errors between filterbands and the errors 
due to the instrumental noise when computing the tempera- 
ture errors and the optical depth errors in the next sections. 



3 PHYSICAL PROPERTIES OF THE WARM 
AND COLD COMPONENTS 

3.1 Far-infrared emission 

3.1.1 Striped region 

We have determined the colour variations along the stripe- 
like region, in the illuminated and shadowed sides (see Fig. 
bottom). For both sides, the correlations Iu{X) — 1^(200) 
shown in Fig. |1| present a linear behave within the scat- 
ter at wavelengths of 120 and 150 ^m. This is described as 
unimodal behave in Paper I. At shorter wavelengths it is 
observed that each correlation is bimodal for the shadowed 
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Table 3. Slopes in the regressions of the striped region observed 
using ISOPHOT. Uncertainties in the slopes correspond to sta- 
tistical variations. 



Region 


^ref 

(pm) 


Slopei 


Slope2 


Shadowed 


60 


0.10±0.02 


0.00±0.02 




70 


0.19±0.01 


O.lOitO.06 




100 


0.40±0.02 


0.3±0.1 




120 


0.55±0.01 






150 


0.910±0.007 






200 


1 




Illuminated 


60 


0.118±0.009 






70 


0.255±0.009 






100 


0.48±0.01 






120 


0.629±0.007 






150 


0.89±0.01 






200 


1 





region: the colour variation is modelled by two connected 
linear ramps with a wavelength dependent turnover point. 
The observed flattening in the emissions towards higher 200 
/xm emission indicates a decrease in the colour ratio at higher 
column densities. 

The slopes of the correlations for the unimodal and bi- 
modal cases are determined from a least-squares linear fit. 
The regression fits are overplotted on the correlations in Fig. 
|1] and the slopes are listed in Table |3 Slopei and Slope2 
stand for the slope values of the correlations before and af- 
ter the turnover emission point. 

3.1.2 Squared region 

We have determined the emission maps of the warm and 
cold components using the second approach described in 
Sect. 2.4.3. Emission maps of the warm component at 60 
and 100 /im are shown (see Fig.^top-left and top-right, re- 
spectively). Emission maps of the cold component at 100 
and 200 /xm are shown (see Fig. |S| bottom-left and bottom- 
right, respectively). The Pearson correlation coefficients of 
the different correlations are listed in Table 0] The number 
of independent pixels (Nip) used in the computation is also 
tabulated. The PCCs corresponding to the warm and cold 
components individually are generally better than the PCCs 
of the full emissions. 

3.2 Colour temperature 

3.2.1 Average colour temperatures in the illuminated and 
shadowed sides 

We have determined the average colour temperatures in the 
shadowed and illuminated sides of LDN 1780 from the cor- 
responding colour-corrected slope of correlation 1^(150) — 
7,/(200), assuming that the spectral energies distributions 
are well represented by a modified black body with power- 
law emissivity /3 = 2 as performed in Paper L We have 
overplotted the fitted modifled black body on the spectral 
energy distributions (SEDs) of the shadowed region (see 
Fig. |7| left) and illuminated region (see Fig. |7| right). In 
the shadowed region we observed a bimodal behave, with 



a flattening in the correlations /.^(A) — /i/(200) at higher 
200 fj,m emission. Two modified black bodies with a power- 
law emissivity /3 = 2 have been considered: one which fits 
to the /^(150) — 1^(200) correlation; other one that fits to 
the /i,(120) —1^(200) correlation, which is representative of 
the denser region (using slope2 at 120 /im of Table |^. The 
colour temperatures of the fits are 16.1 and 15.5 K, respec- 
tively. It is observed that the fltted models for the shadowed 
region reproduce within the errors the entire SED for both 
datasets, and therefore no significant excess is observed at 
any wavelength apart from 60 /im in the low density region. 
An uncertainty in the temperature of 0.2 K has been deter- 
mined from the statistical and responsivity variations. 

In the illuminated region a colour temperature of 
15.9±0.2 K is obtained from the colour ratio derive from 
150 and 200 /im emissions. It is observed an emission excess 
at all wavelengths including 100 and 120 /im with respect 
the fitted modified black body with T=15.9 and /3 = 2. 

3.2.2 Colour temperature of the cold component 

A colour temperature map of the cold component (see Fig. 
ISJ was determined from the 7^(100) and 7^(200) maps that 
were obtained using the second approach described in Sect. 
2.4.3. For each independent pixel, the emissions were colour- 
corrected and fitted to a modified black body with (3 = 2. 
The temperature map might present a small gradient with 
temperatures increasing from Ril5.8 to 17.3 K along the NW- 
SE axis for the region enclosed by the 9 MJy sr~^ isophote 
of 7^(200) (see Fig.|H]l, outside which the temperature errors 
significantly increase. The mean temperature is 16.3 K. 

The error in temperature is determined from the in- 
strumental errors of the IRIS 100 Atm and ISOPHOT 200 
/^m filter-bands, as well as the relative errors and the zero- 
level uncertainty (see Sect. 2.4.4). The errors in tempera- 
ture are ±0.5 K in the area closed by the isophote 7^(200) 
above 9 MJy sr~^. The temperature errors at lower emis- 
sions increase to few K and we have not included them in 
our analysis. 

3.2.3 Colour temperature of the warm component 

The colour temperature map of the warm component was 
determined from the 7™ (60) and 7™ (100) maps. The carriers 
in the wavelength range between 60 and 100 /im are gener- 
ally VSGs and "classical" BGs. However, I™(QQ) was built 
up to be free from VSGs and cold BGs (see Sect. 2.4.3). For 
each independent pixel, emissions were colour-corrected and 
fitted to a modified black body with f3 = 2. The mean colour 
temperature is 25 K. We determined errors in temperature 
of ±1 K with the same criteria applied for the cold compo- 
nent but using the errors in the 60 and 100 /im emissions. 
When assuming a lower value of /3, temperature increases. 
For example, for l3—l, the colour temperature of the warm 
component would be 30±1.5K. 

3.3 Optical depths of the warm and cold 
component 

The optical depth of the cold component at 200 nm 
(see Fig. El left) has been obtained using the relation 
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Table 4. Pixel-pixel correlations of the far-infrared emission before component separation (all), and for the warm and cold component. 
PCC refers to the Pearson correlation coefficient and Nip corresponds to the number of independent pixels. Between brackets the 
results for the region with 72(200) 9 MJy sr^-*^. 
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Figure 6. Top-left: Contours in MJy sr ~^ corresponding to the 60 /im emission from the warm component; top-right: 100 ^m emis- 
sion from the warm component; bottom-left: 100 iim emission from the cold component; bottom-right: 200 ^m emission from the cold 
component. 



TS{200){x,y)^I^{200)(x,y)/B^{T), where the Planck func- 
tion BviT) is determined at 200 /im, using the colour 
temperature T derived in Sect. 3.2.2. The map presents 
a peak at 7.3to'g x 10~ , whose location is in agreement 
with the maxima of 1,5(200) and Av- The mean errors are 
Ar^(200) 



The optical depth of the warm component at 100 /im, 
T^{1QQ) (see Fig.|^ right), has been similarly obtained from 
/J(IOO) and the Planck function at 100 /^m and tempera- 
tures derived in the previous section. We have computed 
an error of Ar™(100) =tfoi- The optical depth at 100 ^.m 
can be converted into optical depth at 200 ^m by dividing 
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Figure 7. Left: SED in the shadowed region of LDN 1780. Two modified blaek bodies are overplotted: one with the colour temperature 
fitted to the /^(ISO) — /^(200) correlation; the other one fitted to the /,y(120) — /^(200) representative of the denser region; Right: SED 
in the illuminated region. 



C(IOO) by (200/100)^ assuming /3 = 2. We note that the 
optical depth of the warm component is few tens lower than 
the optical depth of the cold component. 

3.4 The extinction-infrared relation: an 

estimation of the emissivity of the warm 
component 

The correlation 1^(200)- Ay (PCC=0.76) presents a higher 
scatter for low emission with /^(200) below ~ 10 MJy 
sr~^. At this low emission level there are some satellite Ay 
peaks surrounding the Ay maximum that correlate well with 
r™(100), i.e., the infrared tracer of the column density of the 
warm component (see figures |21 and |^ right). This indicates 
that Ay also gauges the column density of the warm compo- 
nent in the external regions. Therefore the total extinction 
can be written as the sum the warm and cold components 
(see also Cambresy et al. 2001): 

Ay^A^y + AZ (5) 

The extinction of the warm component Ay can be re- 
lated to r™(100) using the next expression: 

AX^^r^imT-^^^ (6) 
Assuming that the ratio ^'^ 4"'°"* is constant, it is possi- 

V 

ble to derive Ay from Eq. 6 and subtract it from Eq. 5 to 
determine Ay- 

The value of iL^^"* has been obtained maximising 

V 

the correlation between and the residual Ay=Ay — Ay. 
We obtained ~1 10~* mag^^, which corresponds 

V 

to ^"j^oo) ^0.25 10"" mag-^ for (3=2. This value implies a 



range of A™ between 0.01 and 1.3 mag. For 13=1, it is de- 
rived that = 0.5 10"" mag~^ Fig. ITOl fng/trt shows 

the 7^(200) - Air correlation with PCC=0.80. It is found 
that PCC=0.85 if considering only values of 1^(200) above 
9 MJy sr"^ 

An averaged ratio Ii{2QQ)/A''y = 12.1±0.7 MJy sr"^ 
mag~^ is derived from the 1^(200) — Ay correlation using 
the bisector method (Isobe et al. 1990). The slope is the 
same when using Ay instead of Ay. It is interesting to note 
that the zero level calibration established in Sect. 2.4.4 is in 
agreement with the origin of the 1^(200) — Ay correlation. 
Note that if adding an offset of 0.3 mag to Ay (see Fig. 1101 
right), the 1^(200) — Ay will pass though the origin (i.e., 
7;;(200)(A^ = 0)=0). This offset of 0.3 mag is consistent 
with the scatter in the relation 1^(200) — Ay and also with 
the uncertainty in Ay (see Sect. 2.2). 

From the correlation t^(200) — we have derived that 
r^(200)/A^ = (2.0 ± 0.2)10"'' mag"\ which is similar to 
the value in the DISM, but few times lower than the values 
found in other translucent regions (see Paper I). 

3.5 Relations of the warm and cold components 
with the HI excess emission, W12 and Wvi 

We have observed that the warm component is located in 
the illuminated side of LDN 1780 as well as the HI excess 
emission (Mattila & SandeU 1979). 

The monoxide line integrated W12 and Wi^ do not cor- 
relate with the warm component. The correlations with the 
cold component are very good (see Fig. II H and Tablej^J. The 
correlation between VK13 and t^(200) and also with 1^(200) 
are especially good with PCC=0.96. The morphologies of 
r^(200) and Wu are very similar (see Figs. Oandl^l left). 

The ordinary least squares bisector method has been 
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Ay (mog) A\ (mog) 

Figure 10. Left: 7^(200) versus Ay; Right: 7^(200) versus Ay. Solid lines show the corresponding least-squares fitting, and dotted lines 
indicate the corresponding upper and lower 3a limits. 



used to determine t^{200)/Wi3 from the correlation 
r^(200) - Wi3. We obtained that t^(200)/Wi3 = (5.2 ± 
0.2) X 10"* K"^ Km"^ s. An offset of Ar^(200) = (0.91 ± 
0.07) X 10"" is obtained. 

3.6 Discussion 

3.6.1 Separation of the warm and cold components 

LFHMIC95 already pointed out that the profile (along the 
direction of the Upper Scorpius OB association) of 25 fj,m 
emission in LDN 1780 shows a different distribution with 



respect to the 12 and 60 fim omissions. In the classical pic- 
ture, this could be interpreted as due to a stronger emis- 
sivity of the VSGs, which mostly emit at 60 fim. We have 
alternatively interpreted that the spectra energy distribu- 
tion in LDN 1780 is mostly due to the existence of cold 
BGs together with warm BGs and VSGs illuminated by an 
enhanced local interstellar radiation field. The warm and 
cold components have high temperatures compared to other 
moderate density regions (see Paper I) and consequently 
the bulk of their emissions is displayed towards shorter 
wavelengths. In the illuminated side of the stripe, it is ob- 
served emission excess at 100 and 120 fim over a modified 
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Table 5. Correlations of W\3 and W12 with /i/(200), /^(200) and t^(200). PCC refers to the Pearson correlation coefficient and Nip 
corresponds to the number of independent pixels. W13 and W12 are positive values. Between brackets the results for the region with 
7= (200) ^ 9 MJy sr"!. 
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Figure 8. Colour temperature iiiaj) and contours (in steps of 
0.5 K) of the cold component. Pixel size of 4' X 4'. The uncer- 
tainties ±0.5 K. The white thick contour marks the 9 MJy sr"'^ 
isophote of 7= (200). 



black body with the colour temperature derived from the 
/i/(150) — 7i/(200) correlation. This excess at such long wave- 
lengths can not be due to VSGs and supports the existence 
of an additional BG component. 

In previous studies, the 60 /im emission has been used 
as template to derive the excess emission at 100 /im (Lau- 
reijs, Claxk & Prusti, 1991; Abergel et al. 1994) and also 
200 /xm (see Paper II). The correlation A7^(100) - 7i,(200) 
(PCC=0.70) derived from our first approach is not as good 
as the /;;(100) -1^(200) correlation (PCC=0.98) obtained 
from the second approach (see Sect. 2.4.3). Different mean 
colour temperatures are derived, with ~13K and ~16K for 
the first and second approaches, respectively. The low tem- 
perature of ~13K is likely due to the fact that 7^(60) con- 
tains not only emission from VSGs but also from warm BGs. 
The second approach to separate the warm and cold compo- 
nents is intended to sort out this problem. The S/N of the 
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Figure 11. Correlation W13 — r^(200). Mean errors arc indicated. 
A total of 34 independent pixels with size 6' x 6' is shown. 



emission at 25 iim. is high enough to be used as template of 
the VSG component. 

The goodness of the second approach to disentangle the 
different dust components is supported by the good agree- 
ment between the temperatures obtained for the squared 
field map of the cold component and the shadowed side of 
the stripe-like region. It is also supported by the very good 
correlation between W13 and t^(200). When using the first 
approach the optical depth of the cold component is not so 
well correlated with W13. The A7^(100) - 7^(200) correla- 
tion is not as good as 7^;(100) — Iv{2QQ). In addition, the 
warm component presents a distribution which is similar to 
the HI excess emission, facing the illuminated side of the 
cloud 

Regarding the properties of the warm component (see 
Sects. 3.2.3 and 3.3) it is difficult to derive its emissivity. We 
have found a ratio ti/'" (IQQ) / Ay which is significantly low. 
The temperature could be also somehow overestimated since 
it could be that Iv{2QQ) contains some contribution of the 
warm component, which is not considered in Eq. (4). How- 
ever, the mentioned agreement between the temperature do- 
rived for the stripe in the shadowed region and temperature 
in the corresponding region in the squared field supports the 
consistency of our results. 
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3.6.2 Properties of the cold component 

The cold component exhibits a temperature of «16 K, which 
is higher than expected at the UV-shielded densities of 10^- 
10^ cm~^. This requires either an extra heating rate through 
the cloud or a change in the properties of the grains with 
respect those of the DISM. 

The average value of 7^(200)/Ay = 12.1±0.7 MJy sr"^ 
mag^^ of the cold component in LDN 1780 is below the value 
of the diffuse interstellar medium with 7,^(200)/A\^=21.1 
MJy sr~^ mag~^ and T=17.5K, and also below the values 
of the moderate density regions of the sample in Paper I. 
These regions are generally colder than LDN 1780 and have 
a cold component with an enhanced far-infrared emissivity. 

The average ratio r^(200)/4^ = (2.0±0.2) x 10"* 
mag~^ is slightly lower than the value of the DISM 
(r^(200)/4v=2.4xl0"'' mag~^) and the same than in 
the warmer (T=18.9±0.7 K) G90. 7+38.0 (2.0±0.8)xl0"* 
mag~^) (see Paper I). When considering values of indepen- 
dent pixels, we find a variation in t^{200)/Ay from 1 10~* 
mag~^ to 3 10 ~* mag"^ for temperatures of 17.3 and 15.8 
K, respectively. A small enhancement in the emissivity of 
the cold component is therefore observed. This could be 
related to grain coagulation as observed in other translucent 
clouds (see Paper I). 

3.6.3 Properties of the warm component 

The temperature of the warm component of 25±1 K (/3 = 2) 
is higher than the average value of the DISM (17.5 K). The 
warm component, which is surrounding the cold component, 
it is certainly affected by the anisotropy of the radiation 
field. The maximum of the far infrared emission of the warm 
component is located in the illumination side of the cloud. 

For (3=2 we obtained r™(200)/A^=0. 25x10"* mag"\ 
This value is significantly lower than the values for the cold 
component as well as the DISM and G90. 7-1-38.0. All these 
regions are significantly colder than the warm component 
in LDN 1780. This supports that the properties of the dust 
grains of the warm component are different from those of the 
DISM as well as the cold component. This could be result 
of the enhanced radiation field in the surroundings of LDN 
1780 with respect the DISM. The observation of the emis- 
sion of all the area LDN 1780 at intermediate wavelengths 
between 100 and 200/im would help to a better extraction 
of the warm component. 

Note that /3 could be different for the warm and cold 
components (see also discussion in Paper II). An unreal- 
istic value of /3=1 for the warm component implies that 
r™(100)/yl^ decreases ~50%. If /3 ~1.5, then T~27K and 
the ratio (100)/ Ay drops ~20% with respect to the value 
for 13=2. A further investigation is out of the scope of the 
data presented here. 



4 lONISATION PROCESS 
4.1 Relation Av ~ /^(Hq) 

The Ha emission in LDN 1780 (see Fig. [HI left) varies be- 
tween ~ 1 and 4 Rayleigh, and its error is ~0.4 Rayleigh 
(Finkbeiner 2003). The extinction Av correlates very well 
with the Hq emission, with a Pearson correlation coefficient 



of 0.71 (see Fig.[ni right). There is a linear relation between 
the Ha emission and Av that remains till Av ~2 mag, and 
then flattens. The correlation is observed in regions with an 
abundant molecular gas content, which is surprising since 
the UV radiation field from the Galactic plane and the Sco- 
Cen OB association cannot come in through so deep (see 
Sect. 4.2.1). 

Such a good correlation between the Ha emission and 
the extinction in LDN 1780 strongly suggests that the 
ionised gas and the dust share a significant emitting volume 
and are quite uniformly mixed within it. The flattening in 
the Ha emission at high Av (~2 mag) is likely result of dust 
extinction (see Fig. 1121 right). It is therefore suitable to cor- 
rect the observed Ha emission 7^,oi)s(Ha) from extinction 
using the following equation: 

r m^^ /^.obs(Ha) r(Ha) 
^'^(H") ^ (l-e-^(H^)) 

The optical depth at the wavelength of Ha is ob- 
tained from our near-infrared extinction map and the 
Cardelli, Clayton & Mathis (1989) extinction law which 
gives Aiia/Av = 0.72 so that r(Ha) = 2.5 log e x 0.72 x Av. 
In the surroundings of the cloud, with r(Ha) 0, the equa- 
tion (6) g 1 VGS 1 1/ (Ha) « 7,y,obs(Ha); toward the peak of the 
cloud /i,(Ha) « r(Ha) 7^,oi,s(Ha). 

The Av — Iv(iia) plot, with the extinction-corrected Ha 
is shown in Fig. 1121 (right). Using a ordinary least squares 
bisector method (Isobe et al. 1990), we obtained a ratio 
7^(Ha)/Av=2.2±0.1 R mag"\ This is twice the value of 
1°'"' (Ha) / Av . Note there is an offset in the relation with 
7^(Ha)=1.4±0.15 R at Ay = mag. 



4.2 Enhanced ionisation rate 

The disk-like distribution of emitting gas in the Galaxy 
produces an emission 7^ (Ha) ~ sin(b) Ra-yleigh (Reynolds 
1984), which corresponds to 1.7 R at the latitude of LDN 
1780. 

The observed relation between 7y(Ha) and Av indi- 
cates that the ionisation and recombination processes are 
even produce in very deep regions of the cloud, with op- 
tical extinctions of ~2 mag. This is likely due to a source 
of ionisation that can penetrate much deeper into the cloud 
than the UV and the X-ray radiations and produces a nearly 
uniform ionisation rate in the cloud. 

The observed ratio I,^{}ia)/Av gives us the possibility 
to estimate the ionisation rate from the following equation: 

I. (He) 
^ ~ N{H+H2) 

which results to be ~10^^^ 7 s"^. We have used 
N(^H + H2)/Av = 1.87 x lO^^ mag"^ cm'^ (Bohlin, Sav- 
age & Drake 1978). Although this corresponds to a very high 
ionisation rate for the intermediate density cloud LDN 1780, 
it worth noting it is only a lower limit obtained by assum- 
ing that every collision with an high energy particle yield 
a Ha photon emission. Other clouds in the LDN 134 like 
LDN 169, LDN 183 and LDN 134 also present Ha emission. 
The origin of the ionisation source is discussed in the next 
sections. 
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4-2.1 Influence of the X-ray and UV radiation fields 

LDN 1780 as well as LDN 134, LDN 183 and LDN 169 are 
likely influenced by an anisotropic radiation field (THLM95, 
LFHMIC95). The southern sides of these clouds face the 
nearby UV sources C, Oph (an 09.5V star at a distance of 140 
pc from the Sun) and the Upper Scorpius OB association. 

LFHMIC95 studied the variation of the dust and molec- 
ular gas tracers along a line parallel to the direction of the 
UV radiation field (assumed toward Upper Scorpius OB as- 
sociation) crossing the centres of LDN 1780 and LDN 134. 
They pointed out some differences between these clouds. For 
LDN 1780: i) no detection of dense and opaque clumps is ob- 
served, although THLM95 found a weak emission of C^*0, 
ii) the emission is significantly brighter at 25, 60 and 100 
lira emission compared to LDN 134 (also LDN 183); iii) the 
60 and 100 /im profiles arc asymmetric (as in LDN 134), 
but similar, and peak at the same position as the carbon 
monoxide lines, while in LDN 134 there is an offset for the 
W(^^CO); iv) /,.(60)/J^(100) and 1^(12)//,. (25) decreases 
toward the shadow side, while in LDN 134 they arc cither 
nearly constant or present a decrease toward higher densi- 
ties; v) the CO temperature peak is highly symmetric, con- 
trary to LDN 134; vi) the lines are narrower on the shadow 
side. 

THLM95 suggested that only the external region of 
LDN 1780 might be affected by the local radiation field. 
They estimated that the radiation field in the illuminated 
side of LDN 1780 could few (2 or 3) times higher than the 
average radiation field at its Galactic latitude because of 
the presence of the OB association. THLM95 concluded the 
asymmetric UV illumination from the Scorpius-Centaurus 
association and from the Galactic plane has infiuenced the 
morphology of LDN 1780, resulting in the offset of the HI 
emission maximum, the blue-shifted (out-flowing) molecu- 
lar gas envelope and the enhancement of the VSG abun- 
dance toward the side facing the Galactic plane. THLM95 



claim that the offset found between the radiation temper- 
ature peaks of CO and ^''CO and the comet- like morphol- 
ogy cannot be explained by the UV asymmetry. They con- 
cluded from the study of the morphologies of CO(J=1-0) 
and ^^CO(J=1-0) that the offset between the maximum of 
the CO temperature peak and the geometrical centre can 
be explained as result of a shock front that passed through 
the cloud. The core of ^^CO and the evaporating lower den- 
sity envelope of gas pointing to NE correspond to the "head" 
and the "tail" of the cometary nebula, respectively (see next 
Sect.). 

The maps of Ha emission and Ay show a good corre- 
lation, and no asymmetry evidence is observed between the 
shadow and illuminated sides. The observe Ha enhancement 
is entirely related to LDN1780. We also observe that Ha 
emission is present in the other clouds of the LDN 134 com- 
plex. However, LDN 1780 is the only cloud of the LDN 134 
complex where ERE emission has been reported (Chlewicki 
& Laureijs 1987). Smith & Witt (2002) have interpreted the 
ERE in LDN 1780 as result of photoluminescence of sili- 
con nanoparticles. They observed a relatively high quantum 
efficiency of the ERE process (~13%) as well as a long wave- 
length emission peak at 7000 A. They argued that given the 
relatively high electron density and low UV photon den- 
sity within the cloud in comparison with the DISM, the 
overall ionisation equilibrium tends more towards neutral- 
ity, and then the silicon nanoparticles are predominantly 
neutral, with an increment in the abundance of the largest 
particles, and then the overall quantum efficiency increases 
and the emission peak shifts to longer wavelengths. We have 
checked the Ha emission maps and far-infrared emission of 
the 23 objects (DISM, LDN1780, reflection nebulae, HII re- 
gions) in the sample of Smith & Witt (2002) finding that 
there are only few of them with a clear counterpart in Ha. 
LDN 1780 is one of the few of them with a correlation be- 
tween Ha and far-infirared emission. We have also noticed 
that other clouds of the LDN 134 complex as the dense LDN 
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183 present Ha in regions with a very low emission at 12 and 
25 ^m. 

We conclude that the asymmetry of the radiation field 
is mainly affecting the warm component, which mostly con- 
tain warm BGs. The warm component is likely associated to 
the HI excess emission and it is surrounding the cold com- 
ponent (see Figure ISJ. The cold component, consisting of 
cold big dust grains, is well shielded from the UV radiation. 
The observation of Ha emission at high column densities 
where CO emission is also observed supports the existence 
of a ionisation flux that can penetrate very deep into the 
cloud reaching the inner regions of it with Ay ~2 mag. The 
UV radiation field from the Galactic plane and the OB as- 
sociation significantly decreases for Av ~0.3-0.5 mag and 
can not explain the observed CO and Ha emissions. The 
X-rays dominates the cosmic ionisation for smaller column 
densities. Therefore, both ionisation sources cannot give ac- 
count of the observed ionisation rate in LDN 1780 and can 
be rejected out. 

4-2.2 Is it possible that the Ha emission in LDN 1780 has 
been produced by a shock? 

THLM95 pointed out the importance of a shock front of 
~10 Km s~^ in the morphology of LDN 1780. Based on the 
idea of de Geus et al. (1989), THLM95 suggested that an ex- 
panding HI shell originated from the Scorpius-Centaurus OB 
association could be eventually fragmented forming clumps 
that were reached by a shock front (of 10-15 km s~^) of 
supernovae originated in the Upper Centaurus-Lupus sub- 
group of Sco-Cen association. As result, some cometary 
shaped clouds as LDN 1780 were formed. 

Assuming that LDN 1780 was formed as result of a 
shock front, is it possible that the emission in Ha has been 
induced by a shock? There is an evolution in the Ha surface 
brightness, /,^(Ha), of a region that undergoes a shock and 
compress to form a molecular cloud. During the initial phase, 
when the gas temperature is still high (~ 1000 K) the typical 
surface brightness of Ha is ~ 0.005 (0.009) R for a shock 
speed of 10 (20) Km s~^. In general, once H2 is formed 
(dominated by self-shielding or dust shielding depending on 
the shock speed, see Bergin et al. 2005) nearly all hydrogen 
is molecular (only about 1% by fraction in atomic form). 
Only a small fraction of H"'" is produced, which will give Ha 
when it recombines. This should drastically reduce the Ha 
emission once the gas is (mostly) molecular (Edwin Bergin 
2005, private communication). Higher values of the shock 
speed front yield higher values of /,^(Ha) at the first stages 
of the shock, but /,^(Ha) always drops by several orders 
once H2 is formed. The timescale of CO molecular cloud 
formation is not established by the H2 formation rate on 
to grains, but by the timescale for accumulating a sufficient 
column density (Ay 0.7 mag; Bergin et al. 2005). 

4-2.3 A plausible solution: confinement by self generated 
MHD waves of cosmic rays 

We interpret the observed ionisation rate for LDN 1780, 
5 ~ 10~^^7 s~^, as due to an enhanced cosmic ray flux. 
The observed value of ^ is around one order of magnitude 
higher than the standard value for the ISM (^ ~ 10~^^ s). 



Our determination of ^ depends on the assumed standard 
ratio Av /Nn^- This ratio could depend on the environment 
conditions; a factor 2 for the uncertainty is probably conser- 
vative. 

Enhanced values of ^ have been suggested for the DISM 
by McCall et al. (2003) from H3" measurements, but it is by 
first time found in a translucent cloud. 

Low-energy cosmic rays of ~100 MeV dominate the ion- 
isation fraction and heating of cool neutral gas in the inter- 
stellar medium, especially in dark UV-shielded molecular 
regions (Goldsmith & Langer 1978). The observed Ha en- 
hancement could be explained as result of a confinement 
by self generated MHD waves of low energy cosmic rays 
(Padoan & Scalo 2005). There are other alternatives that 
require further investigation as multiple magnetic mirrors 
that could lead to cosmic ray variations (Cesarsky & Volk 
1978). 

The study of Padoan & Scalo (2005) indicates the pos- 
sible enhancement of the cosmic ray flux for densities below 
~5xl0^ cm~"^. Although this value depends on different as- 
sumptions (see Padoan & Scalo 2005 for discussion), it is 
in good agreement with the value found for the ^'^CO core 
of LDN 1780, with a density of 6 x 10^ cm"^ according to 
THLM95, and 10^ cm"^ for LFHMIC95. Our relatively high 
temperature (~16K) for the cold component, whose emis- 
sion correlates very well with Av and ^^CO, also supports 
the presence of a source of ionisation like the cosmic rays 
that go much deeper than UV photons and contribute to 
the heating of the cloud. 

The interaction of cosmic rays with molecular clouds 
can produce Ha by several ways. The more direct one would 
be Ci? -f H2 H -I- H+ -I- e" where the H+ emits Ha 
by recombination with an electron. The result of the reac- 
tion also contains a neutral H which can be excited and 
lead to an Ha photon emission. The collision of an en- 
ergetic particle with molecular hydrogen does not neces- 
sarily produce an hydrogen atom. The reaction could be 
CR -I- H2 ^ HJ + e~ which reacts quickly with another H2 
(Lepp 1992): H+ + H2 ^ H+ + H followed by the classical 
dissociative recombination + e^ H-|-H-|-Hor H2+H. 
All the hydrogen atoms resulting from these reactions can 
potentially emit an Ha photon when they are produced in an 
excited state. A fully description of the processes involved in 
the cosmic ray interactions in the ISM to give account of the 
Ha emission has to be carried out to conflrm our interpreta- 
tion. In particular, the efficiency of each reaction mentioned 
above and its probability to produce an Ha photon is un- 
known. The cosmic ray rate we have estimated from the 
Ha surface brightness is believed to be a conservative lower 
limit. 



5 CONCLUSIONS 

The analysis of the data shows the following: 

(i) The Ha emission is well correlated with the extinction 
map, with an evident increase up to optical extinctions of 
~2 mag. There is a deep penetration of the ionisation fiux 
that produces Ha in regions with a molecular content, but 
also a likely infiuence of the ionisation fiux on the balance 
of the chemical abundances in the densest regions. Cosmic 
rays of ~ 100 MeV can reach the inner regions. 
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(ii) The warm and cold components of large dust grains 
in LDN 1780 have been separated. The warm component 
is surrounding the cold component, and it is mainly in the 
illuminated side of the cloud as the HI excess emission. The 
cold component is associated to the ^'^CO core. 

(iii) The warm and cold components have an average 
colour temperature of 25 K and 16 K (assuming (3=2), re- 
spectively. The optical depth of the warm component is few 
tens lower than the value of the cold component. The ob- 
served colour temperatures arc relatively high for moderate 
density regions (see Paper I), which support the presence of a 
large penetration depth ionisation source (i.e, cosmic rays). 
The cold component presents a ratio T2qo/Av that varies 
between 1 10^^ and 3 10^^ mag^^, which arc significantly 
above the values of the warm component and around the 
value of the diffuse interstellar medium (2.5 lO""* mag~^). 

(iv) The warm component consists of warm BGs that are 
influenced by the asymmetric radiation field, with has a ma- 
jor component from the Sco-Cen OB association and the 
Galactic plane. The UV radiation field, with a penetration 
depth much lower than the cosmic rays, might aifect the 
outer regions of the warm component, which might be asso- 
ciated to the HI excess emission {Av up to 1 mag). The inner 
regions of the warm component could be partly associated 
with CO; the blue-shifted (out-flowing) CO component be- 
ing likely affected by the UV radiation field, although ^^CO 
correlates better with the cold component. The cold compo- 
nent consists of cold large dust grains that are shielded from 
the UV radiation field, but not from the incident cosmic ray 
flux. 

(v) The value of the observed ionisation rate ^ ~ 10~^® 
7S~^ for LDN 1780 is relatively high in comparison with 
the values for dense regions (fewx 10"^*^ 76-"^) and for the 
warm (T~8000K) and cold (T>50K) neutral medium (~1.8 
10~^^ ■ys~^), but consistent with the enhancements found by 
McCall et al. (2002) in the diffuse ISM (fewxlO"^^ 7s"^). 

(vi) Recently, Padoan & Scalo (2005) have proposed the 
confinement of low-energy cosmic rays of ~100 MeV by self 
generated MHD waves in the cloud, that can predict en- 
hancements in the ionisation rate of up to 100 times the 
standard value for densities of --^500 cm~^. This value is 
consistent with the densities derived from CO observations 
in LDN 1780. The enhanced cosmic ray flux is likely associ- 
ated to supernovae in the Sco-Cen OB association. 
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